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Abstract 

We explore the phase structure for defect theories in full generality using the gauge/gravity 
correspondence. On the gravity side, the systems are constructed by introducing M (probe) 
D(p + 4 — 2A;)-branes in a background generated by N Dp-branes to obtain a codimension- 
k intersection. The dual gauge theory is a U{N) Super symmetric Yang-Mills theory on a 
{1 + p — A;)-dimensional defect with both adjoint and fundamental degrees of freedom. We 
focus on the phase structure in the chemical potential versus temperature {fi, T) plane. We 
observe the existence of two universality classes for holographic gauge theories, which are 
identified by the order of the phase transition in the interior of the (/x, T)-plane. Specifically, 
all the sensible systems with no defect show a third order phase transition. Gauge theories 
on a defect with {p — l)-spatial directions are instead characterised by a second order phase 
transition. One can therefore state that the order of the phase transition in the interior of 
the (/X, T)-plane is intimately related to the codimensionality of the defect. We also discuss 
the massless hypermultiplet at low temperature, where a thermodynamical instability seems 
to appear for p <2>. Finally, we comment on such an instability. 
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1 Introduction 



Gauge/gravity correspondence [1-4] provides a powerful tool to investigate the dynamics of 
strongly coupled gauge theories. The original formulation [1] conjectures the equivalence 
between supergravity on the "near-horizon" geometry generated by a stack of N coincident 
D3-branes {AdS^ x S^) and the gauge theory (four-dimensional = 4 SU{N) Supersymmet- 
ric Yang-Mills) living on the boundary of AdS^, which describes the brane modes decoupled 
from the bulk. It can be straightforwardly extended to any asymptotically AdS x Ai geome- 
try, M being a compact manifold. This conjectured equivalence is made precise by identifying 
the string partition function with the generating function for the gauge theory correlators, 
with the boundary value of the bulk modes acting as source of the correspondent gauge theory 
operator [3]. 

It can be extended to the case of arbitrary Dp-branes {p ^ 3), for which the world- 
volume gauge theory is again equivalent to the supergravity on the near-horizon background 
generated by the Dp-branes [5]. The dual gauge theory is a (p+ l)-dimensional U{N) super- 
symmetric Yang-Mills theory. Contrarily to the case of the D3-branes, it has a dimensionful 
coupling constant and the effective coupling depends on the energy scale: the gauge theory 
is no longer conformal. More specifically, there exists a frame [6] in which the near-horizon 
geometry induced by the Dp-branes is conformally AdSp-\-2 x 5^"^ [7^9]. In this frame, the 
existence of a generalized conformal symmetry [10] becomes manifest and the radial direc- 
tion (transverse to the boundary) acquires the meaning of energy scale of the dual gauge 
theory [9,11], as in the original AfiS/CFT-correspondence. Moreover, the holographic RG 
flow turns out to be trivial and the theory flows just because of the dimensionality of the 
coupling constant. In the case of the D4-branes, the theory flows to a 6-dimensional fixed 
point at strong coupling: the world-volume theory of D4-branes flows to the world-volume 
theory of M5-branes. 

Gauge/gravity correspondence can be further generalized by inserting extra degrees of 
freedom in the theory. More precisely, one can add a finite number of branes and consider the 
probe approximation, so that the backreaction on the background geometry can be neglected. 
Inserting probe branes introduces a fundamental hypermultiplet in the gauge theory, partially 
or completely breaking the original supersymmetries [12]. 

Here we are mainly interested in the phase structure of the BPS brane intersections at 
finite temperature and finite chemical potential. The phase diagram temperature versus 
chemical potential has been studied in details especially in relation to the D3/D7 system, 
where the probe D7-branes are parallel to the background D3-branes [13-20]. A similar 
analysis has been carried out for the D4/D6 system with the direction of the background 
D4-branes which is not parallel to the probe D6-branes compactified to a circle so that the 
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system has effectively codimension [21]. Considering a system at finite temperature and 
finite ("baryonic") chemical potential means considering a black brane background solution 
and a non-trivial ("electric") profile for the world- volume gauge field on the probe branes: the 
temperature of the system is given by the Hawking temperature of the black brane background 
and the chemical potential is provided by the boundary value of the time component of the 
world-volume gauge potential. 

Along the temperature axis, the system undergoes a first order phase transition (Figure 1), 
which can be seen as a meson dissociation transition [14,17,22]. From the branes perspective. 



Minkowski Embedding Critical Embedding 




(c) Black-Hole-Embedding 

Figure 1: First order phase transition. At low enough temperature (T < T^) the probe 
branes lie outside the black hole (Minkowski embedding). Increasing the temperature, the 
black hole increases in size and the probe branes starts to bend towards it more and more 
until they touch the black hole at T = Xj (Critical embedding). For T > part of the 
brane falls inside the black hole (Black hole embedding). The transition from Minkowski 
embedding to black-hole one is of first order. 

at sufficiently small temperatures, the probe branes lie outside the black-hole background 
(Minkowski embedding - Figure 1(a)). Increasing the temperature, the size of the black hole 
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increases and, consequently, the probe branes start to feel stronger and stronger the black 
hole attraction, bending more and more towards the black hole until they touch the black hole 
horizon at one point (critical embedding) at temperature T = (Figure 1(b)). Continuing 
increasing the temperature part, of the probe branes fall into the black hole horizon (black 
hole embedding - Figure 1(c)). 

Along the chemical potential axis, the system instead undergoes a second order phase 
transition [23]. In the interior of the (/i, T)-plane there has been identified a transition curve 
fi = m (T) : below this curve the system is in a Minkowski phase in which the "quark" density 
is zero, while above it the system is in a black-hole phase with non-vanishing "quark" density. 
This transition curve was numerically found to be first order [16]. However, recently an 
analytic computation by Faulkner and Liu [20] showed that this transition is actually of third 
order. In this picture, in the region fi < m (T) of the phase diagram (Minkowski embedding), 
the DBI-action of the D7-branes provides the dominant contribution to description of the 
brane embedding. For /i > m{T) string worldsheet instantons contribute as well [24], 
creating an instability: the instantons condense and create a neck between the probe D7- 
branes and the black hole so that at /u = m (T) the Minkowski embedding goes over to a 
black-hole embedding (Figure 2). 



Minkowski Embedding Inslanton Condensalion 




(a) Minkowski-Embedding (b) Instanton-Condensation 

Figure 2: Worldsheet instanton effects. In the region fi > m the instanton corrections 
to the DBI-action becomes relevant given their dependence on ~ e '"'^("^"^^ « Their 
contribution becomes dominant so that they condense following a neck until the black hole, 
sending the system to the black hole phase. 

A natural question to ask is how universal are these phase transitions and which type 
of phase transitions these systems allow to investigate. The intersecting brane construc- 
tions provide an arena to construct model which may points towards a deeper theoretical 
understanding of QCD-like [25-29] and condensed-matter-like [30-37] features. 
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In this paper, we investigate the phase structure in the plane chemical potential versus 
temperature for BPS intersecting brane systems, which are constructed by introducing a stack 
of M probe D(p + 4 — 2fc)-branes in a background generated by Dp-branes (M <^ N). 
The parameter A; (A; = 0, 1, 2) indicates the codimensionality of the intersection, i.e. the 
number of the spatial directions of the background branes along which the probe branes do 
not extend. The dual gauge theory is therefore a {p+ l)-dimensional U{N) Supersymmetric 
Yang-Mills theory with a (p + 1 — /c)-dimensional defect on which the fundamental degrees 
of freedom propagate. 

We start with considering the massive fundamental hypermultiplet, which means describ- 
ing the embedding of the probe branes through a coordinate of the transverse space. In order 
to introduce the chemical potential, we turn on an electric ansatz for the world-volume gauge 
field. 

At zero temperature one can identify a brane/anti-brane phase and a "black-hole" crossing 
phase. In the first case, the probe brane can extend from the boundary to a minimum 
distance at which the branes turn and go back hitting the boundary again, forming therefore 
a brane/anti-brane configuration. In the second case, the probe branes can extend down 
to the location of the background branes. The latter phase is thermodynamically favoured, 
the grand-potential being negative. The phase transition between these two configurations 
is of second order as in [23]. The presence of a second order phase transition along the 
chemical potential axis is therefore a universal property of the Y)p/T){p + A — 2k) system. The 
previous analysis applies to systems with k = 0, 1. The case A; = 2 is very different and 
subtle: turning on the embedding mode in the transverse space does no longer correspond to 
a mass deformation, but rather it provides a vacuum-expectation-value for the dual operator. 
However, suitably defining the "physical" chemical potential as well as the "physical" mass, 
one can show that also gauge theories dual to these brane constructions fall in this universality 
class. 

At finite temperature, the systems show a transition line in the chemical potential versus 
temperature plane. We analyse such a phase transition analytically following the approach 
proposed in [20]. The equations of motion for the embedding function and the gauge po- 
tential are too involved to be analytically solved. The idea in [20] is to perturbatively solve 
them considering the "quark" -density (suitably rescaled to be a dimensionless parameter) as 
parameter expansion. Such a perturbative analysis allows to obtain an analytical expression 
for the chemical potential which in turn allows to analytically study the order of the phase 
transition. In the case of systems with no defect, the chemical potential acquires the following 
expression up to first order in the quark density Cf. 

H = m{T)+s, (T) - 5, (T) log + O (c^) , (1.1) 
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i.e. the structure found in [20] for the D3/D7 system generahses to more general Dp/D(p + 4) 
systems. 

In the case of codimension-1 defect systems, there is no logarithmic term at first order. 
This is actually crucial. The presence of the logarithmic term in (1.1) implies that the phase 
transition is a of third order. In the case no logarithmic term is present at first order in Cf, 
the system show a second order phase transition. It is thus possible to state that the order 
of the transition line in the chemical potential versus temperature plane is strictly tied to the 
codimensionality of the system under examination: systems with a codimension-0 defect (i.e. 
no defect) show a third order phase transition, while systems with a codimension-1 defect 
are characterised by a second order phase transition. In the latter case, the whole transition 
line is of second order. For A; = 0, the function S2 has a point at T = T^. < where it 
vanishes: this point represents a tricritical point at which the phase transition becomes of 
second order. 

We also consider the case of massless hypermultiplets in the fundamental representation. 
As mentioned before, such description is provided by fixing the probe branes to wrap the 
maximal S^~'' C S"^"^ with the probe branes wrapping in a (p + 2 — A;)-dimensional subspace 
of the (p + 2)-dimensional non-compact manifold. 

For defect theories, this embedding of the probe branes is described by the scalar mode 
x'' = z (p), with being the direction of the background branes along which the probe branes 
do not extend. When such mode has a non trivial profile, the supersymmetries are broken 
since it provides a vacuum-expectation- value to its dual operator. The supersymmetries 
are instead preserved if the scalar mode z (p) has a trivial profile (i.e. no scalar mode is 
turned on). For codimension-0 {k = 0) systems, the super symmetric embeddings are the only 
possible. In [38] massless N = 2 hypermultiplets in = 4 Supersymmetric Yang Mills were 
studied via probe /^(T— 2A;)-branes Avrctppm^ Sill .AdS^ _ ^ x ^S''^ space. In this framework, the 
specific heat and the zero-sound mode for holographic quantum fluids have been discussed, 
i.e. translationally invariant systems at low temperature and finite chemical potential. The 
zero sound mode turns out to persist to all the value of the hypermultiplet mass for such 
systems [39] and was extensively discussed for D4/D8/D8 systems [40]. In [41] the analysis 
of [38] was extended to T)p/T)q systems, where both massless and massive deformations were 
analysed. In real world, there are two main classes of quantum liquids in (1 + (f)-dimensions, 
with d > 1: the Bose and Fermi fluids. In the first case, the low energy elementary excitations 
are given by superfluids phonons with linear dispersion relation and specific heat c„ ~ T'^. 
Fermi liquids have both bosonic and fermionic quasi-particle excitations, with the latter 
dominating at low temperature and fixing the specific heat to scale directly proportionally to 
the temperature: c„ ~ T (Vd > 1). This Fermi-liquid type of behaviour has been observed 
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in holographic gravity duals constructed with background D4-branes [40,41]. The (1 + 1)- 
dimensional case is indeed peculiar and can never be described through Bose or Fermi liquids. 
First, in (1 + l)-dimensions just collective motion is possible. The reason is straightforward 
to understand. The excitations have one-direction only in which they can propagate and 
therefore they would necessarily scatter with other excitations put them in motion. Such 
systems are thought to be described by the Tomonaga-Luttinger liquids [42-47]. They do not 
have quasi-particle peaks and are characterised by two elementary excitations: plasmons, i.e. 
sound waves as response of charge density to external perturbations with speed dependent on 
the interactions, and spin density waves which propagates independently of the plasmons. A 
first attempt to extensively discuss a holographic realisation of (1 + l)-dimensional quantum 
liquids has been made using a D3/D3-system [37]. At low temperature, the specific heat in 
the dual field theory scales quadratically in the temperature, while for Tomonaga-Luttinger 
liquid the specific heat scales linearly if the system is not interacting and scales as ~ j^^iK) 
for interacting systems (the power a{K) depends on the Luttinger parameter K and acquires 
different forms according to the repulsive or attractive nature of the interaction). 

In this paper we focus again on the massless hypermultiplet. As in [41], we notice that at 
low temperature both the entropy density and the specific heat scale with the temperature 

p-3 

as T^-p, which implies that the systems with p < 3 have divergent entropy density and 
divergent specific heat as the temperature approaches zero. Furthermore, the specific heat 
appears to be negative. For p = 3 the system is characterised by a non-zero entropy density 
at zero temperature and the specific heat scaling with the temperature depends on the codi- 
mensionality of the defect: c„ ~ 2^2(3-fc)^ Finally, for p = 4, the entropy density and the 
specific heat scales linearly in the temperature and, therefore, such a system is characterised 
by a zero entropy density and a zero specific heat at zero temperature. While the thermo- 
dynamical behaviour for p > 3 is physically intuitive, it is not the same for p < 3. First, a 
scaling T~"' (7 > 0) for the entropy density and specific heat seems to violate the third law of 
thermodynamics, which requires the entropy density of a system to reach a minimum^ at zero 
temperature (and, consequently, to have zero specific heat). We study the thermodynamical 
stability of these systems. Since we work in the grand-canonical ensemble, we need to check 
the positive-definiteness of the Hessian matrix of (minus) the grand-potential. The Hessian 

^Strictly speaking the third law of thermodynamics (Nernst theorem) states that at zero temperature the 
entropy vanishes - provided that this limit is taken by keeping the other thermodynamical quantities fixed. 
However, there exist systems which Nernst theorem does not apply to or, more precisely, where its statement 
needs to be modified. For example, the previous formulation does not apply to amorphous solids, which are 
not in equilibrium. They have a finite non-zero entropy at zero temperature (holographic quantum liquids 
constructed with background D3-branes have such a characteristic [38]). Therefore, a more general formulation 
of the third law of thermodynamics is that the entropy reach a minimum in the limit of zero temperature, 
which is the one we stated in the text. 
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matrix turns out to be positive definite for p > 3 and negative definite for p < 3. This 
is a signature of a tliermodynamical instability at low temperature for p < 3. The natural 
question to ask is what is the interpretation of such an instability and how it matches with 
the existence of perfectly well behaved backreacted solutions at zero temperature involving 
D-branes with p < 3, such as the Cherkis-Hashimoto solution for D2/D6 systems [48]. We 
will argue that most likely this is a signature of the breaking down of the probe approxi- 
mation. As we just mentioned, the Hessian matrix of the grand-potential turns out to be 
negative definite with just one of the eigenvalues being negative. One might think to try to 
stabilise the system by turning on also a magnetic component for the world-volume gauge 
field strength. This however will just produce an overall factor dependent on the magnetic 
charge and, as a consequence, it does not stabilise the system. The only tunable parameter 
is the number M of probe branes. Increasing it, one obtains a family of potential curves and 
it is possible to keep increasing the number of branes until the Hessian changes sign. We will 
elaborate more extensively on this point in the Conclusion section. 

The paper is organized as follows. In Section 2 we introduce the codimension-fc defect 
theories by discussing their brane realisation. The codimension-A; defect is created by intro- 
ducing probe D(p-|-4 — 2/c)-brane in the background generated by the Dp-branes in a suitable 
way. In Section 3 we consider the probe brane embedding configuration which introduces 
massive fundamental degrees of freedom in the {p + 1 — A;)-dimensional defect. We explore 
the {fi, T)-plane for such systems. We show that along the chemical potential axis (T = 0) 
there is a second order phase transition for any sensible system of codimension k = 0, 1. We 
briefly discuss the case k = 2. We also investigate the existence of a phase transition in the 
interior of the (/i, r)-plane, showing that the order of the phase transition crucially depends 
on the CO dimensionality of the system. In Section 4 we consider the probe brane embedding 
configuration which introduces massless degrees of freedom. We discuss the behaviour of the 
low temperature density entropy and specific heat. We analyse the stability of the system 
by studying the positive-definiteness properties of the Hessian matrix and we show that the 
systems with p < 3 are thermodynamically unstable at low temperature. Finally, Section 5 
contains conclusion and a summary of the results. 



2 Codimension-/^ Defect Theories 



Consider the background generated by a stack of N black Dp-branes in the near-horizon limit 



-) ' [-h,ir)dt'+d-^^]+(^ 



r„\ — 



dr"^ 
hp {r) 



+ r^dnl_^ 



(2.1) 
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where the constant and the function hp (r) are respectively given by 



J-P 



2^)' P r ( ^ ) g^N {a') ^' = d^g^N {a 



7 — p 



7-p 



7-p 



hjr) = 1 



7-p 

_h 



(2.2) 



with parametrising the position of the black brane horizon The dilaton and the background 
{p + l)-form are respectively 

(7-p)(p-3) 



9. 



) Co...p — g^^ 



7-p 



(2.3) 



The coupling constant ^ym of the dual gauge theory is dimensionful (for p / 3) and defined 

by 

5^ 9A2^r-Ha'f-^ ■ (2.4) 
The temperature of the background is given by the Hawking temperature 

S-P 

T = — = — — . 2.5) 

2tt iirrp \rpj ^ ^ 

Let us redefine the radial coordinate according to the following differential relation 



da 



dr 



so that the background metric (2.1) takes the form 



(2.6) 



7-p 

a \ 2 



11+ 



h2 



dt'^ + dlc"^ 



+ (^) ^ hl~^[da' + aW,_p], (2.7) 



7-p p-3 



where the new radial coordinate a has been rescaled 



<7 > , CTh = —, 



a £ [dh, +oo[ 



(2i 



2 7-p 



and the functions h^{cr) is defined as 



7-p 

■> I ^7-p 



(2i 



The position of the black-brane horizon is now parametrized by and the background 
temperature can be rewritten as 



T 



5-p 

7-P fa,,\ ^ 



9-p 

27-P7rr„ 



(2.10) 
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The dilaton and the background (p + l)-form respectively become 



(7-p)(p-3) 

= g^h^ l^j ' , Co..., = g:X l^j . (2.11) 

For later convenience it is also useful to express the background metric (2.1) in a frame 
in which it is manifestly conformal to an AdSp + 2 x -S** " black-hole space for p ^ 5 (the 
so-called "dual" frame) [9,11]. This can be easily seen by redefining the radial coordinate 
according to 

dcf r^-P 5-p 

= —= — , u„ = (2.12) 

and rewriting the line element (2.1) as 

/ ,\2/(7-p) „ 
c;2 — AT^<P\ 



dsi^ = {Ne^j dsi,„ (2.13) 

so that the line element ds^o describes an AdSp+2 x S^^^ black-hole geometry 

dsl = 5^"^ (^ye^^'^\u'[-hAu)dt' + dl^']+[Kiu)r'^ + uldnl_ \ (2.14) 



with 

,7-p 



h,(^) = 1- = -^r^. (2.15) 

Parametrising the position of the event horizon by u^, the background temperature can be 
conveniently written as 

T = = -^n,. (2.16) 

It is useful to redefine the radial coordinate through the following differential relation 

, , du dp , ^ 

[h,(n)] ^ — = ^, (2.17) 

so that the boundary and the black-hole horizon are now located at p = and p = 1 
respectively and the function hp(ii) can be conveniently written as 

. X h^ (P) .X n7-p 

= Kip) = ^^P (2-18) 

Let us rescale the radial coordinate p as follows 

P - = 4^, (2.19) 

where parametrises the position of the black hole horizon. In these coordinates, the 
background metric, the dilaton and the background {p -\- l)-form becomes 



7-p 



r \ 5-p £^ p-3 f 

dsl =[f] p^Ak'-Hp) 



hi(p) 2 ^ ,-.2' 
hi (p) 



(p-3)(7-p) „7-p 
7- \ 2{5-p) P;:! (p-3)(7-p) 1 / r \ o o^-p 



(2.20) 
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From now on we will work with the coordinates (2.7), unless otherwise specified. 

In the geometry (2.7) let us introduce M parallel probe D(p + 4 — 2/c)-branes (M <C N) 
according to the following intersection configuration 
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intersecting the background Dp-branes along p — k directions {a;''|^_^ and thus forming a 
defect of codimension k (k = 0, 1, 2). The probe branes wrap an internal (3 — fc)-sphere 
g^-k ^ The transverse space can be parametrized as 

dsl = dp^ + {de'^ + sin^ OdOl _ , + cos^ ed^i _ ^ + ,) 

dsl = de'^ + dy^ + Q'^dnl_, + y^dnl_^^^, 

where the two parametrisation in (2.21) are related by y = cjcos^, g = usin^. The above 
configuration for the brane intersections ensures that the systems are BPS at zero temperature 
and therefore no stability issues arise for the ground state. 

The presence of probe branes introduces hypermultiplets in the fundamental representa- 
tion propagating on a (1 + (p — /c)) -dimensional defect. Given the presence of a codimension- 
defect, the embedding of the D(p+4 — 2A:) branes in the Dp-brane background can in principle 
be described through two functions = z[p) (for k ^ Q) and either = 6{p) or y = y (p) 
dependently on the parametrisation (2.21) chosen. Let us comment on these two classes of 
embeddings. 

Fixing the position of the probe branes in the ((1 + (p — A;))-dimensional non-compact 
submanifold at z = 0, the probe branes embedding can be controlled by a scalar mode 
which can be turned on by requiring that a coordinate in the transverse space has a non 
trivial profile. As mentioned earlier, the embeddings can be parametrized either through 
9{p) or y{g), dependently on the parametrisation of the transverse space (2.21). In this class 
of embeddings, the space-time distance between the background branes and the probe ones in 
the transverse space appears as a parameter which is related to the mass of the fundamental 
hypermultiplet. Thus, analysing this class of embeddings is equivalent to consider massive 
hypermultiplets propagating in a (1 -|- (p — A;))-dimensional defect. 

This physical interpretation holds for A; = 0, 1. In the case of codimension-2 intersections, 
the gauge theory on both branes stays dynamical and can be viewed as two gauge-theories 
coupled through bifundamental hypermultiplets living on a (1 + (p — 2))-dimensional defect 
(the D3/D3 system was discussed in [49], while the Higgs branch for Dp/Dp systems was 
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analysed in [50]). Furthermore, the coefficient of the normaUzable mode determines the vev 
of the operator dual to the embedding function, implying that the space-time separation 
between the background Dp-branes and the probe Dp-branes is no longer a parameter but 
rather provides a vev for a dynamical field. 

The second class of embeddings describes the embedding of the (p + 2 — /c)-dimensional 
submanifold in the non-compact (p + 2)-dimensional manifold A4p + 2 keeping the 

position of probe branes in the transverse space fixed to wrap the maximal sphere ''. This 
corresponds to set the mass of the fundamental hypermultiplet to zero. It can be parametrized 
through the coordinate x'' = z{p), with 9 = 0. The embedding mode z is related to the 
vev of its dual operator O,. If it is has a non-trivial profile, the operator O, acquires a 
non-zero vev breaking the supersymmetries. The supersymmetries are instead restored if the 
embedding mode is constant {z = 0). In the case the probe branes fill the whole non-compact 
manifold Mp + 2, obviously it is not possible to turn on such a mode and the description of 
the massless excitations is necessarily supersymmetric. 

For the time being, let us keep both of the two embedding functions and consider an 
electric component for the world-volume gauge field strength F^g 



-f (p) dt A dp. 



(2.22) 



This means that the gauge field on the boundary field theory couples to a [/(I) C U{M) 
current. With such an ansatz, the probe branes are described through the DBI-action which 
acquires the form 



So 



(p- 



- 2k) 



-MT, 



MTo^p + 2)K I d^+"eii_h' 



p+l-2fc 



7-p 



,5-p 



n 5-p 



(2.23) 



The second line of (2.23) has been written by using the solution for the dilaton (2.11) and 
the usual factor (27ra') in front of the gauge field strength has been absorbed in F2. The 
induced metric on the D(p -|- 4 — 2A;)-brane world-volume metric is 



ds^ 

""'p + 6 -2k 



-^dt^ + dx^ 
hi 



+ 



7' 



P-3 
7-p . 



7-p o 5-p 



H r, 



dQ'+['-^ 

a 



''-P p-3 
2 v,7-P^2 



(2.24) 



Q^dn 



2 

3 - fe' 



where x indicates the (spatial) coordinates {x^]^_^ on the defect, the constant A/^. is = 
g~^Yol {S^"*^}, the prime ' indicates the first derivative with respect to the radial coordinate 
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Q and fj^ = + ■ The action (2.23) depends on the embedding function z (p) and the 
gauge field / (p) through their first derivatives only. There is therefore one first integral of 
motion related to each of them 

ll-p-2fc 

V - (2.25) 



o \ z 



Notice that, in the case of a black hole embedding phase for the class of embeddings 2, 
the regularity condition at the horizon, fixes the first integral of motion to be zero and 
therefore the embedding function z must have a trivial profile. In order to have a non-trivial 
profile for the embedding mode, one would need to introduce a magnetic component for the 
world-volume 2-form so that the action for the probe branes has also a Wess-Zumino term. 
The Wess-Zumino action would induce an extra term in the equation of motion (2.25), whose 
value at the horizon fixes c^. We will not discuss the case of the presence of a magnetic 
component for the world-volume gauge field. This means that the regularity condition at the 
horizon forces the system to be super symmetric. 

The first integral of motion is related to the charge density n, which is defined as the 
canonical momentum conjugate to / (p) evaluated at the boundary: 

dL 

n = hm — = MrD(^ + 2)A4(27ra')c^, (2.26) 

Q ^ 00 OJ 

where the factor 27ra' has been restored. For later convenience, let us rescale the coordinates, 
the gauge field / and Cj so that they are dimensionless: 

{t, 1?, Q, y- /} ^ {t, 1?, Q, y- /} , c, ^ Lf-'^c,. (2.27) 

With such a rescaling, the position of the black brane horizon gets parametrized by the 
dimensionless quantity = c^/L*. 

We will discuss the two classes of embeddings separately, i.e. we will discuss both the 
massive and massless degrees of freedom. 

3 Massive hypermultiplet and (/i, T) phase diagram 

Let us now fix the position of the probe D(p-|-4 — 2A;)-branes in the non-compact (1 -|- (p — k))- 
dimensional submanifold {z = 0) and let us consider their embedding in the transverse 
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space, which has been parametrized through the angular coordinate 6 (p). It is actually more 
convenient to parametrise differently the embedding, using a function y (g) of a redefined 
radial coordinate g according to (2.7) and (2.21). Let us write explicitly the induced metric 
on the world-volume of the probe branes 

7-p 



dsi 



h2 



+ 



r X P-3 

a 



l+{yf)dg^ + e'dnl_, , (3.1 



with = g'^ + y"^ and the dilaton given by 



9sh ' 



(7-p)(p^3) 
P-3 / (T \ 4 



a 



7-p 



(3.2) 



where the coordinates are dimensionless. For this class of embeddings, the action of the probe 
branes is given just by the DBI-action 



Su 



(p + 4 - 2k) 



p+l-2fc 



o5-p 



h2 



(3.3) 



As noticed earlier, in the most general case, the action (3.3) depends on the electric gauge 
potential / (g) just through its first derivative, so that there exists a first integral of motion Cf 
related to it. In the limit of zero temperature (i.e. in any point of the chemical potential axis 
in the phase diagram (/x, T)), the action also depends on the embedding function y through 
its first derivative only, so that there is a further integral of motion Cy, beside Cf. In the most 
general case, the equations of motion are 



f'{6) 



1 + (y') 



^/ 5-p 
V,7-p 



cj + g 



3-fch 7~p 



y" 3 - fc / , 2 ^ sy' -y 

l + {y'f Q ^ \<y ) f^^h_h+ 



3 + A;-p + (4-/c) 



a 



7-P 



(3.4) 



^,'-2(3-.)f^)"'4^^ 

In the following subsections, we explore the whole phase diagram for this class of systems, 
including also the simplest cases. 



3.1 The chemical potential axis 

As mentioned earlier, in the limit of zero temperature there is a first integral of motion Cy 
related to the embedding function y, beside Cf. The equations of motion (3.4) have the 
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following simple form 
/' (^?) = 



^2{3-fc) + _ c2 



y' {q) 



g2{3-k) + ^2 



(3.5) 



They can be easily integrated out. Notice that we need to distinguish three cases, according 
to the sign of cj — c^- In this case this difference is negative, the probe branes can extend 
from the boundary to a minimum distance = (c^ — cj)^^'^^^ ^\ where the branes turns 
back and hit the boundary again. The system is therefore a D(p + 4 — 2k) /T){p + 4 — 2k) in 
a Dp-brane background. In this phase, the chemical potential is 



dg 



g2{3-k) + c2 



2(3 - k) _ o 



-B 



J) 



2(3-fe) 



k 1 



2(3- A:)' 2 



(3.6) 



while the embedding function can be explicitly written as 

I c 
yig) = / dg " = 

( 2-k 



B 



V2(3 -ky 2 




k 1 



' 2(3- A;)' 2 



(3.7) 



In (3.6) and (3.7), we set /(a^hj = = y{gmin). Notice that the expressions (3.6) and 
(3.7) for the chemical potential and the embedding functions are valid for k = 0, 1. For the 
codimension-2 systems, the integrals (3.6) and (3.7) are divergent as ^3 — > oo. We will come 
back to this issue later. For the time being, we focus on the case k = 0, 1. 

Taking the limit g ^ oo of (3.7) (in this limit the incomplete Beta- function vanishes), 
one obtains the "quark" mass m in terms of the first integral of motions Cf and Cy. It is easy 
to invert this expression, together with the result for the chemical potential (3.6), to obtain 
the first integral of motions in terms of m and fi: 



2(3 - k) 



B 



1 



2~k 
2(3-fe) ' 2 



3-fc 



fi [m 



2\ 2 



2(3 - k) 



B 



2-k 1 
2(3-fc) ' 2 



3-fc 



m [m 



2\ 2 



The condition cj — < implies necessarily that m > /x for k even. For k 



(3.8) 
1, it is 



straightforward to realize that the expressions (3.8) are meaningful just for in the region 
m > /i (the region m < —fi is not physical). 

In the case cj — > 0, the probe branes can extend to g = 0. Such a configuration 
corresponds to a black-hole crossing phase, in which the chemical potential and the embedding 
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function are expressed by 
dg 



yip) 







^^2(3-/1.) +c2-c2 2(3-A;) (^2_^2)2^ V 2(3 " /c) ' 2(3 - fc) ^ ' 



2 - 



1 



2(3 -A;) (,2_,2 
- B 



(c2 -c2)2{3-fc) 

2 - /c 



B 



2-k 



2(3 -A;)' 2{3-k) 



1 



(3.C 



^2{3-fc)' 2(3- A:)' 2(3- fc). 

At the boundary, the embedding function provides the "quark" mass, which is provided by 
the first term in the second expression in (3.9) (the incomplete Beta-function vanishes). The 
first integrals of motion C/ and Cy can be easily expressed in terms of m and /i: 



2(3 - k) 



B 



1 



2~k 

2(3-fc)' 2(3-fc) 



3-fc 



2(3 - k) 



B 



2-k 1 
2(3-fc) ' 2(3-fc) 



3-fc 



2-k 

m (fx — m ) ^ 



(3.10) 

Similarly to the previous case, it is straightforward to see that equations (3.10) hold in the 
phase diagram region fj, > m. 

Let us now turn to the thermodynamics of our systems in these two phases, by considering 
in particular the grand-potential $7. It is given, up to a sign, by the renormalized on-shell 



action 



where 



D(p + 4 - 2k) 



lim 



(p + 4- 2k) 



D(p + 4 - 2fc) 



(p + 4 - 2k) 



(p + 4 - 2feX 



(p + 4 - 2k) 



Mk / dg 



+ 'S'd(p + 4 - 



„2(3-fe) 



2k) 



(3.11) 



^2(3-fc) + c2 - c2 



(3.12) 



4-k 



(the overall factor Af^ is a redefinition of Mk including J d^^'' and ^ is for cj — > 
or g^in if cj - cl < 0). 

Integrating the expressions (3.12), the grand-potential in the two phases is given by 

(4-fc) 

/„,„„. \ f 

ai (/i^ — 771,2) ■ 



n 



(4-4) 



2(3-fc) 



2(4-fc) 



-B 



2-k 7-2fc 
2(3-fc) ' 2-fc 



(„2 ^2\ 2(3-fc) / 
IS-^/j ^ ' R / 2-fc 1 

2{3-fc)(4-fc) -° \^2(3-fc) ' 2 



4-fc 

2 



a2 (m2-^2^^' 



/J, > m 



fj, < m 



(3.13) 



^The holographic renormahzation of probe D-branes in non-conformal background was extensively discussed 
in [51]. 
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where ai and a.2 are two positive constants. 

Some comments are now in order. Notice that the phase ( "black- hole" crossing phase) 
for which cj — Cy > is thermodynamically favoured with respect to the phase ("brane/anti- 
brane" -phase) for which this difference is negative, given that in the first case the grand- 
potential is negative while in the second case it is positive. Let us analyze the derivatives of 
the grand-potential 0: 



an 



ii(4 - k)fi (//^ - m^) 



2-k 
2 



2-k 
2 



OC Cf 



(3.14) 



OC 



a2(4-A;) 



{^i-k)^?-■m? 



oo 



The point fi = m represents a second order phase transition. This is a universal feature 
for gauge theories whose holographic bulk dual is constructed in terms of flavoured branes. 
In this section we showed how this holds for gauge theories with no defect (k = 0) and on 
a defect with (p — 1) spatial directions {k = 1). The presence of this second order phase 
transition is however common to gauge theories on a defect with {p — 2) spatial directions as 
well, which is dual to the other possible BPS brane construction (A; = 2). In this last case 
one needs to take into account some subtleties which we will discuss in the next subsection. 



3.1.1 Codimension-2 systems 

As pointed out, the analysis in the previous subsection holds for codimension-0 and codimension- 
1 systems, i.e. Dp/D(p + 4) and Dp/D(p + 2) intersections respectively. The Dp/Dp systems 
are quite different from the higher dimensional defect theories. The first difference lies in the 
divergence structure of the on-shell action'^. The case we are now analysing is the only one 
in which the embedding mode saturates the Breitenlhoner-Freedman bound which implies 
the presence of a logarithmic divergence [51,52]. Furthermore, the role of the normalizable 
and non-normalizable modes are exchanged, with the coefficient of the normalizable mode 
determining the vacuum-expectation- value of the operator dual to the embedding function, so 
that the brane separation appears as a vev rather than as a parameter. A similar discussion 



■^For the discussion of the D3/D3 system see section 6 of [52]. The generalisation to any Dp/Dp system 
(p < 5) was discussed in [51] 
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applies for the gauge potential f{p). For both y{p) and f{p), the asymptotic expansion near 
the boundary p ^ oo shows a logarithmic term 

y{p) ''=°° m + Cy log p, f (p) p + Cf log p. (3.15) 

The divergences of the action (3.12) are cured by the counter-terms found in [51]. In principle, 
the integrals of the equations of motion (3.5), which define y{p) and f{p), remain divergent 
as near the boundary. One can define p and m respectively as 

^ = - — , m = —. (3.16) 

Through (3.16), one can rewrite the grand-potential in the form (3.13) (with k = 2). There- 
fore, codimension-2 systems show a second order phase transition on the chemical potential 
axis at p = m as well. 

Thus one can conclude that the existence of a second order phase transition on the 
chemical potential axis is a universal feature of gauge theories with a gravitational bulk dual 
constructed via Dp/D{p + 4 — 2k) systems. 

Before turning on the temperature, let us discuss the behaviour of the chemical potential 
and of the embedding function in the black-hole crossing phase in the limit of small density 
Cf. This will turn out to be useful for the finite temperature case, where the analytic analysis 
is performed in such a limit. 

3.1.2 Small density expansion 

Consider the expression for the embedding function (3.7) taking the limit g — > oo, so that 
it provides an expression of the "quark" -mass in terms of the two first integral of motion Cf 
and Cy. Such an expression can be inverted in order to express Cy in terms of the density Cf. 
In the limit C/ — > 0, one obtains: 



2 2 H 



(3.17) 



where k is constant which contains the "quark" mass as well as other numerical coefficients. 

This suggests that a perturbative expansion around Cf becomes subtle in the region of 
the radial axis for which q'^^^~^^ is of the same order of the leading term in (3.17), i.e. 



1 

g c]-". (3.18) 

Therefore, this region of the radial axis can be conveniently studied by redefining the radial 

1 

coordinate as r = c?"*" g. Keeping this in mind, it is possible to integrate the equation 
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for the gauge field (3.5) to obtain the chemical potential in the small density limit at zero 
temperature 



o 3 — k 

1 + — ^c^-" + O ( c^-" 



(3.19) 



Given equation (3.19), one can expect that the expansion parameter for the chemical potential 



2 

(as well as the embedding function) is c^"*" rather than other powers of the "quark" -density 



Cf. 



3.2 Chemical potential-temperature plane 

In this section we consider the systems at finite temperature and with chemical potential. 
For convenience we rewrite the equations of motions here 



f{Q) 



-■f 5-p 
hi- 



+ 



1 + {y'Y 



3-k I j^^i ^"^ ~ y 



3 + k-p+{4-k) ( — 
a 



7-p 



(3.20) 



^,'-2(3-fc) M'^^ 



In the black-hole crossing phase, the regularity condition at the horizon imposes the following 
boundary conditions for the embedding function 



y' {Qh) = — witho-^ = qI+vI, 



(3.21) 



which prevents the last two terms in (3.20) to blow-up as the horizon is approached. From 
the first equation in (3.20), the chemical potential can be written as 



roc ^_ 



1 + (y'r 



(3.22) 



;2 + ^2{3-fc)h^v-p 



In order to study the transition line in the (//, T)-plane analytically, we closely follow the 
small density approach of [20]. 



3.2.1 Small density expansion of the embedding function 

As noted in [20], the small density expansion for the embedding function and the chemical 
potential is subtle. This is due to the fact that the last term in the embedding function 
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equation (3.20) for Cj — > can be considered as an actual perturbation of the equation at 
zero condensate as long as ^^(^"'^^h.,.'^"'' does not become of order O {cj). For such values of 
the radial coordinates, this term is no longer small. Let us divide the radial direction in two 
regions: 

1. Q e]QA, oo[, 

where is a cut-off distance until which a uniform perturbative expansion can be 
considered; 



1 1 



2. r G]Th, A[, 

where r is defined as^ = Cf^''T, and Th and A respectively as = c]^''Ti^ and 

Qa = C;K. 

The idea is to found a (perturbative) solution in both regions 1 and 2 and then matching 
them in the limit Cj ^ 0, — ^ and A — > oo, keeping r finite. Notice that the above 
splitting of the radial coordinate axis does not hold for fc = 2. i Let us start with finding 
the solution for the embedding function equation (3.20) in the region 1, where, for small Cj, 
the solution can be expanded as follows 

oo 

vie) = ^c>,(£>)- (3.23) 

1=0 

The zero-th order term (g) satisfied the equation (3.20) at Cj = with boundary condition 
yo (0) = 1. In a small g expansion, the 0-th order (in Cy ) solution is given by 

2/0 (^^) = 1 + y^'V' + y^'V' + O (g^) , (3.24) 
where the coefficients yo^\ y[^^ and yf ^ are given by 
(2) _ (3 + A:-p) + (4-fc)^iI-P .7_p 



(4-fc) 



-7-p 



(4) _ ^JlJ y^b(p.fc);^^(7"P) 



4(6-A:)(4-A:)3(^l-<Th^^~^^) ^=o 



(3.25) 



9 

„(6) _ "SP ^(P,fe)^«{7-P) 

24(8 - A:) (6 - A:)(4 - kf (\ - a^^^-^^j i=o 

where the coefficients b and c-''''"* are constants dependent on the spatial dimensions of 
the background branes p and the codimensionality of the defect k and are explicitly given 
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in Appendix A. At the next order in Cf, the solution in a neighbourhood of ~ can be 
generally written as 

(S) = ^^zjr + flo + bo log Q + a,,. ,,Q^-'^ + O (^j2(2-fc)^ . (3.26) 

The coefficient of the leading order o_(2_fc) and the zero order coefficient Oq are constant of 
integration which will be fixed by matching the solutions of region 1 and 2. For /c = only 
the coefficients related to even powers of g are non-vanishing, while for k = 1 both even 
and odd powers are admitted. A crucial difference between the cases /c = and A; = 1 is 
the presence of logarithmic terms in the /c = expansion, which do not instead appear for 
k = 1. The coefficient of the logarithmic term bo for A; = is fixed in terms of the integration 
constant a_2 

.7„p 4(6 -p){5-p)- jp^ + 60 - 29)^.!-^ + 4(7 - p)(5 - p)al^'~^^ + 4^,^^^-^^ 

Oo = ; —o «-2- 

2(l-^,^(^-))' 

(3.27) 

In the expansion for k = 1, the coefficient Oi is determined in terms of the integration constant 

'[n-2p) + 6al-P + {ll-2p)af~^^^a_,. (3.^ 



(7-p)-tI-^ 



Let us now consider the region 2, by recasting the equation (3.20) in terms of the inde- 
pendent variable r. For the solution one can consider the following perturbative expansion 



yir) = ^4z,(t). 
Also, the position of the horizon admits a perturbative expansion for small Cj 

oo 

The 0-th order equation is 



i=0 



= ^ + (3 - A:) 



l+(^ 



+ 2 



7-p 



7-p 



7-p 



1_ (£k 



2{7-p) 



Q - p + 



7-p 



(3.29) 



(3.30) 



(3.31) 
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where the dot ' indicates the derivative with respect to r. The conditions at the horizon and 
at the boundary are 

Oh 

To ' (3.32) 



Zn 



horizon: Zo{to) = cj^, 
boundary: (tq) = 1 

In order to match the solutions in region 1 and 2, it is suitable to look in the region 2 for a 
solution at large r 



(2 - k)i 



i=l 



{2-k)i ' 



(3.33) 



where the very first coefficient ("(2 _ fe) is 

C(2 -k) = — 



(2 -k)[l + al^) 

At first order, the perturbative solution shows a singular term, which is of order O (r): 



(3.34) 



^1 {r) 



^" + 00 + tf,, log r + 



-(2 - fc) 



r 



2-k 



+ 



1 



r2(2-fc) 



(3.35) 



where Oq is an integration constant and the coefficients O2 - fe and VOq are fixed in terms of the 
zero-th order solution to be: 



(3+fc_p)+(4-fc)a;^-P ^7 



for /e = 



(3.36) 



(4-fc)(i-af-); 

y 0, for k = l. 

We are now in condition to compare the solutions in the two different regions of the radial 
coordinate. First, it is convenient to rewrite the expansion in the region 2 just in terms of the 
radial coordinate q as well as the density Cj. We will rewrite the expansions both in region 
1 and 2 so that the comparison between the two results becomes straightforward 

a_ 



-(2 - fc) 
^2-A: 



y(/,)U, = l + y(^)^2^0(/)+c, 

C,2-k 



+ ao + bo log + _ fe) Q^~^ + O ( ^>2(2-fc) 



+ 



+ tto + tt)o log Q 



^logc, + o(,2-'= 



(3.37) 

where the coefficient t)2(2 - fe) comes from the equations at order c^, which we have not written 
down explicitly. The matching at zero-th order in Cj fixes some of the integration constants 
and provides some non-trivial consistency check. Notice that this matching provides a non- 
trivial check on the expansion: 

«2-fc|k=o = 

f2-fcL=i = 0, 



(3.38) 
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where the coefficients m (3.38) are provided in (3.25) and (3.36) and identically satisfy the 
equalities in (3.38). The matching at first order in Cf fixes some integration constants: 

o_(2_fc) = C2-k, tt)o = bo, Oo = tto + ^^logcy, (3.39) 

where the second equality becomes the identity = for k = 1. The constant of integration 
Oo is determined by the boundary condition Ui {g)\g ^ ^ ^ 0- 

This small density analysis of the embedding function showed that the appearance of 
logarithmic terms is common to all the systems with no defect while it disappears when a 
codimension-1 defect is introduced. 



3.2.2 Small density expansion of chemical potential 

We can now use the previous analysis to explicitly compute the chemical potential in the 
small density limit. The chemical potential is given by the integral in (3.22) with extreme of 
integration [q/^, oo[ and [^h, Qa] respectively in region 1 and region 2. 

In region 1, the leading term of the chemical potential in a Cj perturbative expansion is 
of order one 



/^Lcg 1 = c/ / dg {q^ + yl) 
J eh 



ll-p-2fc 



Using the perturbative solution (3.24) for y,), it is easy to compute the first orders of the 
chemical potential in a neighbourhood of 



(1) 



(2 - fc) 



Ircg 1 



2-fc 



+ TTXo log £Ia + nXi^lA + {Qa 



(3.41) 



where 



1 



(2 - k) 



ll-p-2fc ' 
7-p 



(4-^)^(l-a5"^))( 



1 + ^l. 



(7-p)V- 



{p,fe)^j(7-p) 



j=0 



and the last term fC {qa_) in (3.41) is finite and defined by 



(3.42) 



8A 



y = vo 



(2 - k) 



■- fA 



rno log £)a - mi£»A 



(3.43) 
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In region 2, the small density expansion shows both a zero-th and first order terms 



rcg 2 



dT 



(0)1 



I rcg 2 



(3.44) 



where the zero-th order term, which is fixed by the solution (3.33) for Zq, can be conveniently 
written as 



("^l dT^^ (3.45) 



M = m 

I rcg 2 



2 5 — p 

A ^0 / 7-p , -7-p\ 7-p 



Using the perturbative expansion (3.35) as r — > oo, the zero-th order chemical potential can 
be obtained as an expansion in A^^ 



^ Ireg 2 



1 - a 



7-p 



, ii-p-2fc ' ^ V A2(2-fc) / • 

(2-fc)(l + ar^) ^ 



(3.46) 



At first order in c^, the chemical potential receives contributions from the zero-th order 
solution Zq in (3.33) as well as the first order one Zx in (3.35). In full generality, it can be 
written as 



(1) 



I rcg 2 



dT { 



TO 



7-p ^ 7-p 



7-p ^ 7-p 
Zo - CTh 



7-p I -7-pN 7-p 



+ 



2z,f iu Uo + <5-h 



2zl'al-' + (6 - p)zl~^^) - (z',^'"'^ - a, 



AfA7"P) ;^2(7-p)W 7-p , .7-p\ — 



Zo + C^h 



(3.47) 



The first term in (3.47) can be easily integrated. Notice that in the case of codimension-1 
systems, it is the only term contributing to /^'^'|^^g 2- Using the expansions (3.33) and (3.35), 
as well as the relations (3.38) and (3.39), in the case of /c = 1 it acquires the form 



(1) 



1 - a'^-P 



I rcg 2 



5-1, "^O- 



(3.48) 



(1 + (7'^-P)7-P 

Therefore one can write down the chemical potential in the small density expansion as follows 



fi = m 



+ i(r)c, + o(c2) 



where s (T) is defined as 



i(T) = lim /C(^J + //'^)| 

£lA-»0 "^<=S2 



(3.49) 
(3.50) 
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which should vanish in the zero temperature hmit. In principle one would need to fix the 
constant do by imposing the boundary condition yi\^_^^ = 0. However, for our purposes, 
i.e. to check the order of the phase transition, it will not be necessary. The function s (T) 
in (3.49) may be both non-zero and zero. In the second case one would need to go to next 
order in Cf , in which case the small density expansion of the chemical potential /i can at most 
acquire the form: 

^ = m + 5,{T)c^j.- Sa (T) log Cf (3.51) 

However, in any case the order of the phase transition will be the same. We will comment 
on this in the next subsection. For the time being, let us focus now on the structure of the 
small density expansion of the chemical potential for system with no defect. 

In the case of /c = we have a non-trivial contribution to Li^^^''^!! from both the terms 

' Ircg 2 

in (3.47). For convenience let us write (3.47) as 

(Dl'^-o = (1)1 + (1)1 (3_52) 

^ Irog 2 Ircg 2 '^D I rog 2 ' ^ ' 

where A^i^'l^^gj indicates the total term in (3.47). The first term in (3.52) is easy to obtain 
and can be written as 

Lcg2 = — T^^h +- -lEF^o- 

4(1 {i+ai-^y-^ 

_ (6 - p)(3 - p) + (27 - 5p)^il-^ + 4ag("-^) .2(7-^) , 

-2(7-p)W, , ,7-p\2f^ " (3-53) 



4 l-<7:^'--M 1 + 



7-p 



4(6 - p)(3 - V) (1 - + (199 - 26p - 
+ ^ log A, 



where, similarly to (3.48), we would not need to fix the coefficient do explicitly. The second 
term in (3.52) instead acquires the form 

(1)1 _ 3-p + 4aiI"P 2 

Irog 2 ~ s 2^ 



2(p^ - 9p + 12) + (53 + 14p - 3p^)aI~P - 2(12 + 5p - /)a^(^-P) _ 4^^^^"^) 

25-3p 

4fl + <TM 



6 - P + (7 - p)<Th ^ + <Th 



,7-p , ;~2(7-p) 



25-3p 
l—p\ 7-p 



al~'P logA + W 

(3.54) 
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where, similarly to /C, W, is defined as 

W = lim + ^^^^^^al~'A'+c(al~A log A, (3.55) 

4(l + al-)'- ^ ^ 

where C (j^h~^^ indicates the term of (3.54) in curl bracket. Notice that the first (divergent) 
term in (3.54) exactly cancels the quadratic divergence (as A oo) in (3.53). Summing the 
contributions (3.41), (3.46) ,(3.53) and (3.54), as well as using the third relation (3.39), the 
chemical potential in the small density expansion can be expressed as 

fi = m + 5,{T)cf- S2 (T) Cf log Cf, (3.56) 

where Si (T) and S2 (T) are given respectively by 

. m r ^ 1 - ^I"" , (6 -p){3-p) + (27 - 5p)al-^ + 4al^'~^^ ^jr-p) ^ ^, 
Si (T) = /C H —ao -grj fTi, + W 

(l + ^I-)^- 4(l-a,?(^-))(l + ^I-)'-^ 

S2(r) = i(x-y-bo), 

(3.57) 

where x and y are the coefficients of log A in (3.53) and (3.54) respectively, and bg is given 
in (3.27). The expressions (3.49) and (3.56) can be written as 

2 2 

^ = m + 5,{T) c^-" - S2 (T) c]-" log Cf (3.58) 

In the next subsection, we will discuss in detail the order the phase transition in the interior 
of the plane (//, T) using the results of this subsection. 



3.3 Phase transitions in the {fi, T)-plane 

In the previous section we were able to obtain a very general analytic expression for the 
chemical potential in the small density limit. Let us emphasise again an important point. The 
charge density n, which is given in (2.26), is proportional to the density Cj. The Minkowski 
and black-hole embeddings are characterised by zero and non-zero n respectively. Therefore, 
working in the small density limit with black-hole embedding boundary conditions means 
focusing on a neighbourhood of the transition curve in the (/u, T)-plane, so that the limit 
— > send us on the transition curve. 

Let us now look in more detail at this phase transition. Since we are in the grand- 
canonical ensemble, the right thermodynamical potential to analyse is the grand potential 
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and its derivatives. In the Minkowski embedding phase {i.e. fj, < m), all the derivatives of 
the grand-potential turn out to be zero: 



an 



~ n 



0, 



0, 



Vs > 1. 



(3.59) 



In the black-hole phase, {i.e. fj, > m), we can use the expression (3.58) for the chemical 
potential (if s(T) = 0), then take the limit C/ — > to go on the transition curve ^ = m and 
match the result with (3.59): 
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k = 



k = l. 



(3.60) 



Already from (3.60) one can infer that, while for k = there is no discontinuity in the second 
derivative of the chemical potential, for = 1 it diverges. For codimension-1 systems, the 
phase transition in the (/x, T)-plane is of second order. 

Let us now consider the third-derivative of the grand-potential 17 for k = 0: 



(T) 



00. 



(3.61) 



, > ™, . = (Si (T) - s, (T) - 5, (T) log CfY 

As a consequence of such a divergence, the phase transition for codimension-0 systems is of 
third order. 

In the case s(T) 7^ for codimension-1 defects, then the second derivative of with 
respect the chemical potential fi acquires the form 



+ 0{c,] 



Cf-»0 



1 



5i (T) ' 



1. 



(3.62) 



Even with the form (3.49), the second derivative shows a discontinuity. 

One comment is now in order. The results (3.60), (3.61) and (3.62) imply that the order 
of the phase transition in the {fi, r)-plane is tied to the codimensionality of the defect theory. 
We therefore identified two universality classes of theories determined by their codimension- 
ality, or equivalently, by the order of the phase transition in the (/i, T)-plane. 

The existence of a zero in the function S2 (T) would implicate the presence for A; = 
systems of a critical point {fi^, T^) at which the phase transition becomes of second order. 



4 Massless hypermultiplet and quantum fluids from non-conformal branes 
backgrounds 

In order to consider massless degrees of freedom on a (1+p— A;)-defect, we need to fix the probe 
branes to wrap the maximal sphere S'^"'' . The probe branes therefore wrap _[_ 2 — fc ^ ^ 
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Figure 3: Phase diagram (/i, T). The axises are defined as fi/m and T/T^. The phase 
diagram shows a first order phase transition at T = along the temperature axis (systems 
at zero baryonic chemical potential), while at zero temperature and finite density all the 
Dp/D(p + 4 — 2k) systems are characterised by a second order phase transition at /x = mo. 
The curve fi = m (T) actually represents a transition curve in the interior of the (/i, T) 
plane. The region fi < m (T) corresponds to Minkowski-like embeddings with n = 0, while 
the region fj, > m (T) corresponds to a black hole phase, with the transition between the 
two phases is of second order for theories with a codimension-1 defect and of third order for 
theories with no defect. 
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We will work directly in the black hole embedding: since the probe branes wrap the black hole 
geometry Aip + 2- k, the induced metric contains a black hole as well. This class of embedding 
can therefore be described through the linear coordinate x'' = z (r)^ for defect theories. If 
the embedding mode z(r) has a non-trivial profile, its dual operator acquires a non-zero 
vacuum-expectation-value breaking the supersymmetries. If the embedding mode z{r) is a 
constant (namely z = 0), the probe branes have a fixed position in the (p+2 — A;)-dimensional 
non-compact manifold, the operator has a zero vev and the supersymmetries are not 
broken. This is the only possible configurations for theories with no defect (i.e. Dp/D{p-\- 4) 
systems). As pointed out in Section 2, regularity condition at the horizon forces the embedding 
function to have a trivial profile. We therefore need to consider the supersymmetric case. A 
similar analysis was carried out in [41]. We will begin with following [38,41] by computing 
thermodynamical quantities such as the entropy density and the specific heat. 

Let us turn on a non-trivial profile for the world-volume gauge field using the ansatz (2.22) 
so that the dual gauge theory has a chemical potential. The DBI action and the equation of 
motion for the gauge potential / (r) are 

5D(p + 4-2fe) = -MTD(p + 4-2fe) j dt (P - '^x dr d' - V e-<^ V-det {^as + F^^} = 

= -MTD(, + 4-2fe)57'Vol{5='-'=} j dt(F-''xdrr'-'^^Jl-[r {r)f (4,1) 

fir) 



^c) + r2(3-fc) 



The contribution to the grand-potential from the fundamental degrees of freedom is given, 
up to a sign, by the renormalized on-shell action 

= - hm [51_,,,+ 51J (4.2) 

(the holographic renormalization of probe branes in non-conformal backgrounds was exten- 
sively studied in [51]). We can directly compute Ofu^ in the low-temperature limit (i.e. 
Th — > 0), which is the regime we are interested in 

/•r-h ^2(3-fc) 

Td(p + 4 _ 2fe) 



■2(3-fc) + c2 2(4 - A;) 



fun 



2(3 -k) + l 



2{3-fc)+l 



l + O 



r 



2{3-k) ' 



h 



+ 



+ -M rj3(p_,_ 4 _ 2fc)A/'p _ fe— - — 



'2(4 - A;)' 



(4.3) 



Here we indicate z as function of r since we will consider the background line element (2.1) 
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where 



(T = 0) 
fun 



is the grand-potential at zero temperature, which is exphcitly provided in 



Section 3.1 (setting Cy 



0) and we rewrite here for future convenience: 



fun 



(4 - A:)a3-'= 



4- A; 



(4.4) 



Furthermore, the constant J\fp _ ^ is a redefinition of A/^. by including the volume of the defect: 
Afp-k = -^kVp-k = -^k J df^^x. Notice that full grand-potential for the system is given 
by the contribution (4.3) from the fundamental degrees and the one from the adjoint ones: 
^adj = —{5 — p)N vjX^-pT 5-p3_ fpj^g latter contribution, as well as the last term in (4.3) are 
independent of the baryon density. Therefore, in order to study the features of the "quantum 
liquids" we need to focus just on the first two terms of (4.3) which will be indicated as (if^n- 
Similarly, we can compute the chemical potential in the low-temperature limit by inte- 
grating the equation of motion (4.1) for the gauge potential /(r): 



^r2{3-fc) + cj 



Mo 



„2(3-fc) 



2 [2(3 - A:) + 1] cj 



+ 



4(3-fe)' 



(4.5) 



with fXg = aCf * being the chemical potential at zero temperature. Notice that the parameter 



1 

expansion in (4.3) and (4.5) is the dimensionless ratio r-i-./cj^'' or, equivalently in terms of 

2 

the temperature and the chemical potential at zero temperature, T"'-p /^Uq. 
The entropy density can be now computed as a function of T and Cf. 
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(4.6) 



The leading term in the entropy density in s (T, c^) scales with the temperature as s ~ T^-p , 
which is independent of the codimensionality of the defect. Notice that actually this term is 



^The quantity is just a constant dependent only on p: 



(7_p)3{7-p) 
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of order one in the conformal case p = 3, while for p < 3 it decreases with the temperature 

3-p 

as s ~ T 5-p and for p = 4 it increases hnearly as s ~ T. This would imply that only 
in the conformal case, the entropy at zero temperature is non-zero. Not only. In the limit 
T — > the entropy density seems to blow up for p < 3. We will comment later on this. 

From the entropy density (4.6), it is straightforward to compute the specific heat Cv at 
constant volume and density (in the limit of low temperature) . 
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(4.7) 



For p = 3 the first term in (4.7) vanishes so that the leading term for the specific heat Cy 
scales with the temperature as r2{3-fc) ^g^g [38]). For p = 4, Cv ~ T independently of the 
dimension of the defect in which the fundamental massless degrees of freedom propagates. 
This is actually the same behaviour of Fermi liquids. It is interesting to notice that the 
leading term for the specific heat Cy becomes negative for p < 3. This would be indeed a 
signature of a thermodynamical instability in the canonical ensemble. However, our analysis 
is carried out in the grand canonical ensemble where the thermodynamical stability of the 
system is guaranteed if and only if the Hessian of (minus) the grand potential is positive 
definite 
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(4.8) 
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where the elements {rirr, il^p, f^T^) ^^mmI Hessian (4.8) are exphcitly given by 
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where M = MTj^(^p + 4_2k)J^p-k- In order to have the stabihty of the system guaranteed, 
the Hessian TC (tlf^^ in (4.8) needs to be (semi)-definite positive, which means that the 
eigenvalues need to be positive (or at most zero). This is also equivalent to requiring the 
non- negativity of all the principal minors of Tl ( Clf^,^ I . Being a 2 x 2 matrix, there are just two 



principal minors: and the determinant \7{\ of the Hessian. From the explicit expression 
in (4.9), the smallest minor f^yj-- becomes negative for p < 3. As far as the determinant 
is concerned, it has the following form: 



\n\ 



5 — p 



7 — p 



5-p Cf _ 



P- 



3 ^9 + p 



Ak 



5 — p 2(5 — p) cj 



+ 



(4.10) 



As for the minor f^TX; the determinant (4.10) has the leading term in the low temperature 
expansion which becomes negative for p < 3. Therefore, the Hessian turns out to be negative 
definite for p < 3. In the conformal case p = 3, the first term in (4.10) vanishes and the 
new dominant contribution is always positive. Finally, for p = 4 all the terms in (4.10) 
are positive. The Hessian is therefore positive definite for p > 3. The positive-definiteness 
of the Hessian for the cases p = 3, A insures stability. Moreover, the entropy density at 
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zero temperature is either finite {p = 3) or zero p = 4, in agreement with the third law of 
thermodynamics. The specific heat turns out to behave as ~ T'^(^~^^ (p = 3) or ~ T {p = 4). 

A comment is now in order. The statement of the negative-definiteness of the Hessian 
is a local statement. This means that, in principle, the systems could tend to a stable con- 
figuration. Moreover, the brane configurations analysed in this section show a very peculiar 
behaviour for entropy density (4.6) and specific heat (4.7) at low temperature for p < 3. For 
convenience, let us rewrite it below: 

2 _2_ 3-p 

~ 2 / Att \ 5-p _3;iP ~, n — 3 / 47r \ ^-p T ^-p 

.,(T.c,)=M—[—) „r .-,+..., ,,=AA2^|-^(— j — +... 

(4.11) 

The behaviour (4.11) implies that both entropy density and specific heat increase as the 
temperature approaches to zero. This seems to violate the third law of thermodynamics 
according to which the density entropy reaches a minimum value as the temperature ap- 
proaches to zero. As a consequence, the specific heat should vanish in the same limit. This, 
together with the negative-definiteness of the Hessian (which, at the end of the day, is a 
consequence of the behaviour (4.11)), implies that the configuration analysed here are not 
thermodynamically stable for p <,3. However, we need to recall that zero temperature back- 
reacted solutions do exist for p < 3 which are well-behaved, like the D2/D6 solution found 
in [48]. It is therefore natural to ask how the appearance of this low-temperature instability 
connects with the existence of well-behaved backreacted solutions. One explanation to this 
question can be provided by the fact that these systems can become stable once one goes 
beyond the probe approximation, i.e. it is the the backreaction which stabilises them. The 
only tunable parameter available is the number M of probe branes. One can thus think to 
increase the number of T){p -|- 4 — 2fc)-branes until the Hessian becomes positive definite and 
thus the systems stabilise. This would imply a modification of the potential curve and there- 
fore the effect of increasing M can't be seen from our computation, given that the potential 
curves of the probe case and of the backreacted one are different. In other words, tuning M 
one gets a family of potentials until stability is reached at the backreaction, when the probe 
approximation breaks down. 

One might argue that we are just computing thermodynamical quantities which refer just 
to the fiuid (i.e. to the fundamental degrees of freedom), and that our stability analysis 
does not take into account all the degrees of freedom. This can be done by starting from 
the full grand potential, which contains both a contribution from the adjoint degrees of 
freedom and from the fundamental ones (including the density independent term in (4.3)). 
The contributions from O^dj scales with the temperature as T -'-p , while the last term in 

2 4-fc 

(4.3) as T^~. It easy to see that such a contribution will not affect the leading behaviour 
at low temperature of any of the relevant thermodynamical quantities, such as the density 
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entropy, the specific heat and the Hessian of the grand potential. Thus the instabihty we are 
observing is not indeed due to not keeping into account the whole grand potential. 

In principle, looking at the eigenvalues of the Hessian, one might think that the system 
can be driven to a high temperature stable point. In this case, the contribution from the 
adjoint degrees of freedom becomes more and more important as the temperature increases. 
Simultaneously the probe branes get heat up and can acquire enough stress-energy to even- 
tually backreact. The probe approximation may break down anyway. It seems reasonable 
to think that the backreaction can stabilise the system, and the negative-definite Hessian, 
together with an apparent violation of the third law of thermodynamics, is a signature of the 
breaking down of the probe approximation. 

We would like to stress that a deeper analysis of such an instability is indeed needed, 
since our arguments do not provide a robust proof of the nature of such an instability as the 
breaking-down of the probe approximation. 

5 Conclusion 

In this paper, we investigate the phase structure of Dp/D(p + 4 — 2A;)-systems, where the 
D(p + 4 — 2k) branes are considered in the probe approximation. These systems are BPS 
and introduce a {p + 1 — fc)-dimensional defect in the {p+ l)-dimensional U{N) gauge theory. 
We consider both the probe brane configurations which introduce a massive fundamental 
hypermultiplet and the ones which introduce massless excitations in the fc)-dimensional 
defect. 

The probe D{p + 4 — 2A;)-branes wrap a Aip + 2-k x S^~'' subspace of the background 
geometry A4p + 2 x S^~p generated by a stack of Dp-branes. The embedding of the probe branes 
can be controlled by a scalar mode, which is provided by one of the angular coordinates in 
the transverse space. The separation between the stack of probe branes and the background 
ones is a parameter in the theory related to the quark mass. In this setting we consider the 
system to be at finite temperature {i.e. non-extremal black hole background) and at finite 
chemical potential, by turning on a non-trivial profile for the gauge field on the probe branes 
world- volume. On the chemical potential axis (T = 0), the system undergoes a second order 
phase transition for A; = 0, 1. In the region fj, < ttt-q the probe branes cannot reach the 
location of the probe branes, but they have an extreme point where they turn back and hit 
the boundary again: the probe branes are in a D{p + 4 — 2k)/T){p + 4 — 2k) configuration. 
For > 7T7-0 the probe brane can reach the location of the background branes, in a black-hole 
crossing phase. The second phase is thermodynamically favoured since the grand-potential 
is negative, while in the brane/anti-brane phase it is positive. The case k = 2 is a bit 
special given its different physical interpretation: turning on a non-trivial embedding mode 
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in the transverse space does no longer correspond to a mass-deformation, i.e. the separation 
between the probe branes and the background one can no longer be interpreted as a mass, 
but it provides a non-zero vacuum-expectation-value for its dual operator. 

In the interior of the {fi, r)-plane, it is possible to identify a transition-curve fj, = m (T) 
between Minkowski embeddings and black-hole embeddings. The order of the phase transition 
across this curve is strongly tied to the co dimensionality k of the defect. In the case of k = 0, 
the result of [20] can be extended to any sensible Dp/D(p -|- 4)-systems: the transition across 
the curve fj, = m (T) is a third order phase transitions. Third-order phase transitions are 
not very common in nature. However, there are some meaningful examples to keep into 
consideration. One example is indeed provided by the Gross- Witten model [53] in which 
a third-order phase transition in the large-iV lattice gauge theory in two-dimensions was 
observed and then extended to four-dimensions. Indeed, the physics of the phase transition 
in the Gross- Witten model is very different from the one of the class of theories discussed in 
this paper. It is a weak-to-strong coupling phase transition which occurs at fixed 't Hooft 
coupling (A = A^) and at large-A^. Here instead, the system is studied at both large-A^ 
and infinite 't Hooft coupling and the transition occurs in the strongly-coupled regime (for 
p >,0). The gauge theories dual to these brane construction define a new class of theories 
with a third order phase transition. 

For such systems, the transition line fi = m (T) represents a third order phase transition 
until it reaches the temperature T = T,, at which the phase transition becomes of second 
order. 

In the case k = 1, the phase transition along fi = m (T) is of second order. This is a direct 
consequence of the fact that, for such systems, the small "quark" -density expansion of the 
chemical potential does not show any logarithmic behaviour at first order in "quark" -density. 

Theories with a codimension-1 defect can be reduced to effective codimension-0 theories 
by compactifying the direction of the background on which the probe branes do not extend. 
An example of such systems was studied in [17,21,26]. One can wonder if and how the 
reduction of theories with a codimension-1 defect to effective theories with no defect affects 
the order of the phase transition along the curve fi = m{T). Most likely, it is the case because 
of the introduction of a new scale provided by the compactification radius. However, we did 
not check this explicitly and we leave it for future work. Indeed, if the compactification has 
the effect of changing the order of the phase transition along fi = m{T), it would strengthen 
our statement of the existence of two universality classes of theories with baryonic chemical 
potential, which are identified by the codimensionality of the defect, i.e. by the order of the 
phase transition. 

The massless degrees of freedom on the {p+1 — A;)-dimensional defect are studied by fixing 
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the probe branes to wrap the maximal (3 — fc)-sphere in the transverse space. Fixing also 
the position of the probe-branes in the non-compact directions provides a supersymmetric 
description of the system. In principle one can turn on a non-trivial profile for x'' = z (p) 
(indeed this description is valid for A; / 0). In this case, the operator would acquire a 
non-zero vev and the supersymmetries get broken. We consider the supersymmetric case and 
focus on the low temperature properties of such finite density systems. Following [38,41], we 
compute the density entropy and the specific heat in such a regime. Both those quantities 
turn out to scale with the temperature as 



It is interesting to notice that the powers in (5.1) become negative for p < 3. Indeed, this 
type of behaviour seems to be at least counter-intuitive. For the third law of thermodynamics 
one would expect the specific heat to vanish as the temperature goes to zero and the density 
entropy to reach zero as well or, anyway, a finite value. Instead, for the class of systems 
identified by p < 3 both these quantities appear to blow up in the zero temperature limit. 
In order to investigate this issue, we explicitly analyse the stability of these systems. A 
necessary and sufficient condition for the (local) stability of these systems is that the Hessian 
matrix is positive definite. In the case of interest {p < 3), we show that the Hessian is 
actually negative definite, which identifies an instability at low temperature explaining the 
observed behaviour of entropy and specific heat. However, it arises new questions. First, 
what is the nature of such an instability? How can these systems be stabilised? How does 
this result connects with the existence of a zero-temperature backreacted solution such as the 
D2/D6 of [48]? We pointed out that the only tunable parameter in the theory is the number 
M of probe branes. So one can think to increase it until the potential changes in such a 
way to make the Hessian positive definite, at which point most likely the Af D(p + 4 — 2k)- 
branes (p < 3) will start backreacting stabilising the system and the probe approximation 
breaks down. Another possibilities is that the analysis of the eigenvalues of the Hessian may 
suggest that the system can be lead to a stable point at high temperature. In this case, the 
contribution from the adjoint degrees of freedom as well as from the density independent 
term from the fundamental ones will start to be more and more relevant. Simultaneously, 
the probe branes get heat up to the point that they can start to backreact. This would lead 
again to the breaking-down of the probe approximation. It is therefore reasonable to think 
that the backreaction is needed to stabilise the system, and consider it as best candidate for 
explaining the appearance of such an instability at low temperature. Once again, we want to 
reiterate the idea that a more detailed analysis about this stability issue is needed, but our 
observation stresses the necessity of an extensive discussion about the validity of the probe 
approximation as well as the need of a deeper understanding of backreacted solutions and 
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the related physics. 
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A Coefficients for the small density expansions 

In this section we explicitly write down the numerical coefficients in the small density expan- 
sion of section 3.2, where we emphasised as in [20] that such an expansion is subtle and it is 
necessary to split the radial axis in two regions. In region 1, at zero order in C/ the embedding 
function is conveniently expanded in a neighbourhood of g ^ and such expansion turns out 
to involve even powers only of the radial coordinate (3.24). The coefficients of this expansion 
have been written in (3.25) in terms of other coefficients b-*'"''^ and c^'''*'^ which depend on the 
spatial dimensions p on the background branes and on the co dimensionality k of the defect. 
Here we list the coefficients b^'''''' and c^^'''^: 

• Coefficients b-''''"' in yg"* 

b(P, fc) = _(4 _ kf (/ _ (12 + k)p + 9{k + 3)) , 

fc) = 2(4 - kf (/ - 2{k + 8)p2 + {2k^ + 18k + 85)p - 2(9^;^ - 2k + 109)) , 
^(P,k) ^ _2(4 _ (4^3 _ 3(^3 _ ^2k^ _^ 52A; - 52)/ 

+ {3k'^ - 2k^ - 2b2k^ + 1512 - 2068)p - (3 k){2bk^ - 296/c^ + 1224A; - 1564)) , 
fe) = 2(4 - yfc) ((4 - k)p^ + 2(A;2 + k- U)p'^ - {2k^ - 2k'^ + 49A: - l2A)p+ 
-F4(3/c^ - 19A;2 + 95A; - 137)) , 
= _(4 _ kf ((4 - k)p^ + {k^ -Ak + 2A)p - 3(12 - A:)(3 + k)) , 
^(P,k) ^ 4{6-k){4-kf. 

(A.l) 
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• Coefficients c in ^o'^^ 

^(P,k) ^ (6- A;)(4-A:)5(ll-p)(9-p)(3 + A;-p), 

^(P,k) ^ 2{A-kf{9-p) [(16 - + (6A;2 + 17/c - 264)p2_ 

-{k^ + 92k'^ - 347k - 960)p + (15A;^ + 2A2k'^ - 1753A: + 24)] , 
^(P,k) ^ 4(4_/.)3^3_^^_^) ^2(5-A;)p^ + 4(/c2 + 29/c- 256)p^+ 

+ (4A;3 - 156/c2 + 33k + 1258)i92 - (74^^ - 1884A;2 + 7565k - 3970)p- 

-(4 - A:)(340A;2 _ ^g^^j. _^ 9573)] ^ 
cCp,*) = _2(4 _ A:)2 [88p5 _^ ^^^^^^s _ ;^3g^2 _^ 208/c - 2616)/- 

- 2(11^"^ + 18A:^ - 1564A;2 + 3552/c - 17176)/- 

- {3k^ - 714/ + 6150A:^ + 4896/.'^ - 37152/c + 209968)/+ 

+ 2{2k^ - 3147/ + 37726A:^ - 107556/ + 1458884A; + 175420)p+ 
+(203A:^ + 14898/ - 226797/ + 941720A:2 - 1807856/c + 755560)] , 
^(P, fc) ^ 2(3 + A; - / [8(2/ - 35/ + 228/c2 - 648A; + 696)p^- 

- 8(4/c^ - 30/ - 299/ + 3940A;2 - 142898/c + 17504)/ + 

+ (35fc^ - 198/ - 1664/ - 3224/c^ + 197280^^ - 924160/c + 1288320)p^- 

- 4(165/c'^ - 3024A;^ + 21040/ - 89722/ + 372024^^ - 1148640A: + 1484736)p+ 

+ (3069/ - 68754A;^ + 605952/ - 2847248/c^ + 8357440A;2 - 15749056A; + 14671808)] , 
^iP,k) ^ _2(4 _ ^8(4 _ j^-^p5 _ ^55^4 _ gg4^3 _^ 593q;j2 _ i^j^gQj^ _^ 25120)/ + 

+ 2(55A;^ - 424/ - 3884/ + 48176^^ - 173528A: + 237536)/- 

- {73k^ + 422A;^ - 15982/ + 42672A;^ + 4157 Uk^ - 2346864A; + 3647040)p^ + 

+ 2(560A:^ - 7005A;^ + 13600/ + 40732/ + 483088^^ - 3681404/e + 6445680)p- 
-(4187/ - 80082/ + 551079/ - 1930700/c^ + 5556816/c^ - 15048120A; + 20930976)] 
= 4(4_A;)2(3 + A:-/ [2(5-A;)(4-A:)/ + (4- A:)(4/ + ll/e- 306)/+ 
+ (32/ - 496/ + 2921A;2 - 10622A; + 20264)/- 

- (398/ - 7196/c^ + 4654/c2 - 149796A; + 190264)p+ 

+ (1100/ - 22516/ + 162269/ - 504752A; + 629168)] , 
= 2(4 - kf [3(4 - A:) (16 - 3A:)/ + 2(4 - A:) (9/ + 16A; - 516)/+ 
+ (27/ - 206/ + 450/c2 - 8200A; + 33792)^^- 

- 2(136/ - 1896/ + 9704A:2 - 32996A; + 79184)p+ 
+(405/ - 8014/ + 56877A:2 - 182228A; + 322176)] , 

(A.2) 
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+ (131A:2 - 782A; + 2952)] , 
^(p.fe) = 8(4-A;)^(llA;2-24A;-32). 
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